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l. PROJECT  GOALS 

There  has  been  considerable  progress  in  recent  years  in  the 

development  of  quantum  mechanical  methods  for  obtaining  state-to-state 

cross-sections  and  rate  constants  for  simple  chemical  reactions. 

(2  31 

In  addition  to  the  H + H_  reaction,  ’ results  for  F + H„  and  H + D 

(4.5)  Z L 1 

have  appeared.  * These  latter  studies  have  made  use  of  natural 
collision  coordinates  (NCC)^^  and  of  centrifugal  decoupling  approxi- 
mations which  have  been  tested  by  comparison  witli  accurate  close- 

(4  7) 

coupling  calculations  for  H + ’ Such  approximations  make  tractable 

coupled-channel  studies  of  reactions  involving  atoms  other  than  hydrogen. 

The  studies  mentioned  above  have  utilized  an  approximation  to 

the  kinetic  energy  which  is  only  applicable  to  reactions  where  the  most 

(8 ) 

stable  configuration  is  linear.  This  excludes  from  consideration  most 

of  the  bimolecular  reactions  that  are  important  in  mechanisms  of  combustion 
processes.  Such  reactions  may  be  associated  with  radicals  that  are  stable 
in  highly  bent  configurations.  A representative  list  is  given  in  Table  I 
(see  page  32). 

The  primary  goal  of  this  project  is  the  development  of  state- 
of-the-art  computer  codes  capable  of  studying  reactions  involving  non- 
linear intermediates  on  a single  adiabatic  electronic  surface.  Due  to  the 
complexity  of  such  reactions,  stability  of  the  numerical  algorithm  employed 


as 


to  integrate  the  scattering  equations  is  of  prime  importance,  as  weil 
its  efficiency.  Wyatt's  approximation  to  the  NCC  kinetic  energy  will  be 
generalized  to  include  terms  important  for  nonlinear  geometries. 

A secondary  goal  of  this  project  is  the  application  of  the 
scattering  codes  to  processes  of  interest  in  combustion.  Of  particular 
interest  are  the  reactions 


H + 02  ► Oil  + 0 
0 + H2  *■  OH  + H 


which  are  two  important  chain-branching  steps  in  the  - 0^  reaction, 
leading  to  the  formation  of  two  product  radicals  for  each  reagent  radical. 
The  first  reaction  is  highly  endothermic,  and  therefore  slow  at  low 
temperatures.  Knowledge  of  the  temperature  dependence  of  the  rate  constant 
for  this  reaction  is  necessary  for  an  understanding  of  low  temperature 
steady-state  reactions  in  the  - 0^  system,  as  well  as  the  explosion 
regimes . 

In  addition  to  providing  a theoretical  tool  for  calculating  net 
reaction  probabilities,  state-to-state  probabilities  are  obtained  as 
immediate  results  of  such  calculations.  The  role  of  reagent  internal 
energy  in  determining  energy  disposal  for  a reaction  can  therefore  be 
studied  using  the  methods  developed  during  this  project. 
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II. PROJECT  ACCOMPLISHMENTS 

A.  Summary 

During  the  duration  ol'  this  project,  the  natural  collision  co- 
ordinate theory  of  quantum  reactive  scattering  was  generalized  to  reactions 
of  the  form  AB+C  in  which  the  reaction  intermediate  could  be  nonlinear. 

An  approximation  to  the  NCC  kinetic  energy  operator  was  employed  which 
systematically  reduces  to  that  of  Wyatt  when  applied  to  a linear  inter- 
mediate. A scattering  code  developed  previously  was  modified  to  compute 
the  many  additional  terms  required  in  the  kinetic  energy.  The  resulting 
code  is  the  present  state-of-the-art  in  the  application  of  quantum  re- 
active scattering  theory  to  general  bimolecular  reactions. 

The  new  code  was  tested  by  application  to  the  H + 1^  and 
F + reactions,  and  then  applied  to  the  H + 0,,  combustion  reaction.  As 

will  be  seen  below,  many  new  results  were  obtained  which  should  aid  our 
understanding  of  microscopic  processes,  many  of  which  are  difficult  to 
study  experimentally. 

B.  Details 


The  quantum  mechanical  Hamiltonian  of  Wyatt  has  been  generalized 
to  the  case  of  a nonlinear  reaction  intermediate.  We  have  determined  that 
additional  curvature  terms  resulting  from  the  noneuclidean  nature  of  natural 
collision  coordinates  (NCC)  make  important  contributions  to  the  scattering 
matrix  elements  for  large  deviations  from  linearity  in  the  reaction  zone. 

The  role  of  a floating  origin  originally  introduced  by  Elkowitz 
and  Wyatt  lias  been  investigated.  We  have  found  that  a systematic 
choice  of  floating  origin  can  minimize  some  derivative  couplings,  and  makes 
little  difference  in  computed  probabilities.  A choice  of  switching  angle 
can  be  made  which  causes  certain  matrix  elements  to  identically  vanish, 
but  in  reactions  involving  nonlinear  intermediates  can  make  other  matrix 
elements  large. 


The  actual  physical  motion  of  three  interacting  atoms  was  in- 
vestigated by  varying  the  natural  vibrational  and  rotational  coordinates. 

We  have  found  that  a proper  choice  of  NCC  parameters  can  approximately 
uncouple  the  vibrational  and  rotational  motions,  which  simplifies  the 
fitting  of  the  potential  surface.  This  can  lead  to  efficient  ways  of 
choosing  the  points  for  ab-inl t io  calculations  of  potential  energy  sur- 
faces . 

Accurate  fits  were  made  in  NCC  to  several  potential  surfaces, 
including  the  H + 0^  surface  of  Causs.  This  particular  system  is  difficult 
because  of  the  angular  dependence  of  the  reaction  path.  The  use  of  NCC 
makes  possible  the  study  of  correlation  diagrams  relating  vibrational  and 
rotational  states  of  reactants  and  products,  even  when  the  minimum  energy 
path  is  nonlinear. 

Two  versions  of  a scattering  code,  REACTOR,  were  tested,  one  of 
which  utilized  the  R-matrix  propagation  procedure  of  Light  and  Walker, 
the  other  an  initial-value  method  developed  by  Redmon  and  Wyatt.  Both  were 
satisfactory  for  the  H + and  H + 0^  reactions,  but  the  R-matrix  code 
was  found  superior  for  the  F + reaction,  with  its  highly  closed  channels, 
due  to  the  way  in  which  inter-sector  transformations  are  constructed.  In 
all  cases  the  R-matrix  code  was  easier  to  use  because  of  its  inherent 
stability  to  closed  channels.  The  initial-value  code  was  slightly  more 
efficient. 

A study  was  made  of  the  effect  of  variation  in  the  F + bending 
potential  on  the  probability  for  reaction  from  the  ground  vibrational  state 
of  (v=0)  into  the  v=2  state  of  HF.  It  was  found  that  simply  moving  the 
maximum  in  the  rotational  eigenvalue  correlation  diagram  could  alter  the  low 
energy  resonance  found  for  this  process,  and  affect  the  dynamics  of 
population  inversion.  This  illustrates  the  importance  of  reliable  calcu- 
lations of  potential  energy  surfaces  for  nonlinear  geometries. 

Calculations  of  reaction  probability  as  a function  of  initial 
vibrational  state  show  that  reagent  vibrational  energy  is  very  effective 
in  promoting  reaction  in  the  H + 0^  system.  This  is  in  agreement  with  the 
conclusion  of  Gauss,  based  on  classical  studies  of  this  reaction.  An 
important  result  is  that  nonreactive  encounters  of  H with  0^  can  form 
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vibrationally  excited  0 2»  which  will  then  react  readily  with  H atoms. 

Even  at  low  relative  energies,  the  probabilities  for  inelastic  collisional 
excitation  of  0 2 are  comparable  to  the  elastic  scattering  probabilities. 

C.  Publications  Resulti ng  From  This  Project 

The  following  papers  are  in  preparation  for  submission  to 
scientific  journals: 

1.  A Generalization  of  Natural  Collision  Coordinate  Reactive 
Scattering  Theory  to  Nonlinear  Intermediates  (A.  B.  Elkowitz). 

2.  Kinematics  of  Reactive  Scattering  in  Natural  Collision  Co- 
ordinates (M.  J.  Redmon) . 

3.  A Model  Potential  Surface  in  Natural  Collision  Coordinates 

for  the  Reaction  H + 0^  -*•  OH  + 0 (M.  J.  Redmon  and  A.  Gauss,  Jr.). 

4.  Quantum  Mechanical  Investigations  of  the  Sensitivity  of  Re- 
action Probabilities  to  Features  of  the  Potential  Surface: 

The  Bending  Potential  for  F + H2  (M.  J.  Redmon  and  R.  E.  Wyatt). 

5.  Quantum  Reaction  Probabilities  for  the  Reaction 

H + 02  -*■  OH  + 0 (M.  J.  Redmon,  G.  Wolken,  Jr.,  and  A.  B. 

Elkowitz) . 

6.  Dynamics  of  Reactive  Scattering  in  Natural  Collision  Co- 
ordinates (M.  J.  Redmon). 

7.  Quantum  State-to-State  Cross-Sections  for  the  H + 0^  Reaction 
(M.  J.  Redmon). 

8.  Quantum  Dynamics  of  the  F + 1^  Reaction:  Resonance  Models, 
and  Energy  and  Flux  Distributions  in  the  Transition  State 
(S.  Latham,  J.  McNutt,  R.  E.  Wyatt,  and  M.  J.  Redmon,  to 
appear  in  the  Journal  of  Chemical  Physics). 
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111. TECHNICAL  DISCUSSION 

A.  Reactive  Scattering  in  Natural  Collision  Coordinates 

The  coordinate  system  used  .o  describe  the  reaction  dynamics  is 
known  as  natural  collision  coordinates . These  are  a set  of  center  of 
mass,  body  fixed  coordinates  in  which  Euler  angles  0<f>x  are  used  to  orient 
the  three  particle  triangle  in  space,  and  three  coordinates  s,  p,  y are 
used  to  describe  the  shape  of  the  three  particle  triangle.  in  asymptotic 
regions  where  the  separated  atom-molecule  limit  is  valid,  s describes 
translation  of  the  free  atom,  p represents  vibration  of  the  diatomic  mole- 
cule, and  y is  rotation  of  the  molecule  around  the  atom-molecule  axis. 

An  abbreviated  description  of  the  coordinates  as  used  in  the 
current  treatment  is  presented  in  Figures  1 and  2.  Figure  1 shows  only 
those  terms  necessary  to  define  a collinear  geometry  while  Figure  2 indi- 
cates modifications  necessary  to  include  deviation  from  coll inear ity . In 
Figure  1,  translation  is  defined  as  arc  length  s along  the  reference  curve 
(RC),  vibrations,  p,  are  measured  perpendicularly  to  the  RC  (i.e.,  in  a 
plane  of  constant  s),  and  bending  motion,  y,  is  measured  in  a plane  of 
constant  s around  a suitably  chosen  origin  (Figure  2)  such  that  y=0,  it 
denotes  a coliinear  geometry.  As  in  previous  work,  the  RC  is  linear  in 
asymptotic  regions  and  follows  the  arc  of  a circle  with  turning  center 
(TC)  in  the  reaction  zone.  Following  standard  practice,  the  TC  is  located 
so  that  the  so-called  "schizoid  region"  occurs  in  an  area  of  high  potential 
energy.  Asymptotically,  the  RC  is  chosen  to  lie  at  a distance  r = r^  from 
the  appropriate  Z axis.  Unlike  previous  treatments,  however,  no  reaction 
path  (RP  in  1)  is  defined  in  the  collinear  geometry  since  now  the  inter- 
mediate state  may  be  non-linear  (y>0).  Furthermore,  for  systems  of  non- 
linear intermediates,  considerably  more  care  is  needed  in  choosing  the 
floating  origin  (FO)  than  for  systems  with  linear  intermediates. 

The  existence  of  a multi-valued  region  in  the  y=0,  tt  plane  is 
well  known  and  is  solved  by  choice  of  turning  center.  A less  well  known 


trouble  area  exists  in  the  transition  region  tor  i near  90°  where  the  co- 
ordinate system  becomes  overdet  erminod  when  the  switching  angle  a is  not 
equal  to  its  asymptotic  values.  We  have  introduced  a scaled  angle  > de- 
t ined  on  the  intervals  (0 , n/2)  and  (»/2,'>)  such  that  > ranges  only  over 
intervals  (0,  > ) and  (n-  which  eliminates  the  difficulty  (see 

Figure  2). 

l'he  classical  kinetic  energy,  T,  tor  a system  with  a nonlinear 
intermediate  in  this  approximat ion  is  defined  by 


uT;  p‘  + ^ r‘v‘  + ^ R‘x‘  sin  >)xy  + 


Vr‘  sin'x  x " -i±r~  f sin  cosy  a'  x -^r  fx  y + ursin  yA  $x 

2.  y z x 


where 


* - (M  •Mb-Mc/(Ma+Mb+Mc)) 


M , M , are  masses  of  the  three  reacting  atoms  in  any  arrangement 
A D C 


tube . 


T|2  “ (1  + nnl*  + 


dY 

* m ~ Ts 


da 

'l  ds 


n * r - r cos  > 
c 


f - cos  fct-Yl 


m ■ r sin  v 


Ag  - R^Hj  cos  (Q-?)  - (1  + Kn  + h1  n)  sin  yy-a'I  - f 
R - Re2  - rc2  + r2  + 2n  (rc  + sin  (?-Y>). 


3r 


ds 


t he 


, and 
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Therefore,  at  s = 0 t » 0 so  wo  approximate  the  product  l'A„  hy  0 and 

? 

eliminate  the  term  [P  , P ] . Secondly  we  assume  that  in  Q , sin  y -*•  sin 

s Y * 

y^,  where  "0"  denotes  the  "equilibrium"  or  reaction  path  values.  We  thus 
obtain  the  Hamiltonian  operator 


1 

, - 2 


2^  , ,2Q 

P <r  -p)  c 


(rc  - p) 


B 2 Q 


+ 


2n(rc-prQc 


X <rr-  P)QA  + 


o p 


,2„2 


Y l-h  +i-?-  ] — T * Sin  Y ¥ 

2p,  r 2 2 sin  y y Y 

xo 


B sin  y 
2^  r 20  X 


r A 


s 1 


s_  + J-ha.  ~ 


1 Ps  1 


—2  - 2 

BO  B * Q 

o 


(3) 


^ 2 7x2  + — ^y2  + £JZ2  +fBcoty['J  , Jz]+  + 


2p.Q 


2^72  ^ tPY»  sinvl+  + V(s.  p.  y) 

^o 


where  B - p,B 

For  simplicity,  the  potential  energy  may  be  approximated  V(s,P,y)  ■ 

V^(s)  + l^p.s)  + V^(y\s).  In  this  assumption,  is  the  translational 
or  minimum  potential  along  the  non-linear  reaction  path,  is  the  vibrational 
potential  parametric  in  8 and  is  the  rotational  potential.  The  Hamil- 
tonian is  then  partitioned  into  four  parts: 


ms 
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TTF- 


T'r  vtj 
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H + H ..  + H + H 

Crane  vib  rot 


pys 


where 


H 


1 B 2 
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o i VI  -UJ  v/  i i \ - 
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(4) 
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OYS  ? X S 1 1 


PYS  2 a 2„ 
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—2  —2 
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To  obtain  a symmetric  interaction  matrix  the  total  wave  function 

is  rescaled  so  that  Q l/^gl/A 

o 

f = TT^—  *. 


rg 


1/2 


This  leads  to  a Schrodinger  equation  of  the  form 
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where  U( s,  p)  Is  an  effective  potential  with  the  form 
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The  wavefunction  4 is  written  as  an  expansion  in  translation, 
vibration,  and  rotation  wavefunctions  so  that 


jin 
Jo  o o 


ZH«  TJ  l n jtaU)  nH  (MX* s)  »,,<»“> 


where  JM  are  good  quantum  numbers  of  total  angular  momentum  and  Z component 

of  J in  the  body  fixed  system.  Other  indices  are  j.'t.n  which  label  molecular 

angular  momentum,  orbital  angular  momentum,  and  vibrational  energy  respectively 

and  are  good  quantum  numbers  only  in  the  separated  atom-molecule  limit. 

Indices  1 ,1  ,n  label  the  initial  state  of  the  system  before  reaction  begins, 
o o o 

Hn(p;s),  a vibrational  wavefunction,  is  a solution  of  the  equation 

il.  H f {E  -V  (p;  s) ) Hn  - 0 (7) 

2 n 2 n 
Dp  h 

JM 

O':  is  a solution  to  the  rotational  part  of  the  Hamiltonian.  Note  that  in 

J*'  ? - 

the  limit  of  large  atom-molecule  separation,  A “0  and  Q = T|  g so  that  the 

S / \ 


Schrodinger  equation  reduces  to  that  previously  used  by  Wyatt.  As 

the  rotational  wavefunction  0 is  expanded  so  that  in  asymptotic  regions 
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where  + Is  for  products  and  - 1b  for  reactants.  As  the  reaction  proceeds  and 
the  potential  Vj(y;s)  becomes  non  zero  the  rotational  functions  are  expanded  as 


(8) 


The  d’s  are  computed  by  diagonalizing  fi  over  the  rotational  part  of  the 
Hamiltonian  so  that 


Hr°t  wj<tnH 


(9) 


The  rotational  part  of  the  Hamiltonian  in  Eq.  (5)  is  identical 

2 

to  H in  Eq.  (4)  except  it  is  now  multiplied  by  Qq  as  a result  of  the 

. 8 

wavefunction  rescaling.  ^rot  may  be  more  simply  written  as: 
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Equation  (9)  is  multiplied  by  , integrated  over  y,  and 

J ^ 

summed  over  K.  This  leads  to  a set  of  equations  of  the  form: 
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Once  the  terms  in  Equation  (11)  have  been  computed,  this  equation 
may  be  solved  by  diagonalization.  This  yields  both  the  rotational  wavef unctions 

IM  ] 

i4‘, V(d$XY >9)  and  the  energies  W‘  (s). 
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With  SI  and  H having  been  determined,  close  coupling  equations 
J i n 

describing  local  propagation  through  a sector  are 


njt  6JXJ  Vj7  TiotonoJto  (s)  + 

vJ  , <12> 

niJi^i  ^ Tj  tQ^Q  j-tn(s)|  = 0 


where  the  interaction  matrix  is  given  by 
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The  usual  first-derivative  term  does  not  appear,  since  we  are  anticipating 
the  use  of  locally  adiabatic  basis  functions  and  intersector  transformations. 
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Once  the  matrix  elements  V,  (s)  have  been  generated,  the  coupled 

J11  i/yn  B p 

equations  are  solved  via  the  Magnus  method.  Boundary  conditions  are  then 

applied  to  obtain  the  scattering  probabilities. 

The  above  equations  may  be  partially  decoupled  by  the  use  of  so- 

called  J^-conserving  or  coupled-states  approximations.  This  is  an  essential 

procedure  for  applications  to  systems  other  than  H + H,,.  The  calculation 

presented  in  this  report  employed  the  decoupling  approximation  of  Wyatt, ^ ^ 

so  that  calculations  for  J>0  involve  no  more  coupled  equations  than  .1=0. 

Instead  of  the  triangle  inequality,  values  of  j,j  for  J>0  are  determined 

by  the  prescription 


J even  < 


J » odd 


J ■ even 


J odd 


r\  = odd 


.1  “ even 


« = I J-J 


p if  j * J 
[i  if  j > J 


« - J-J i 


p if  J - -1 
'(j  + 1 if  J > 


These  relations  result  in  proper  bifurcation  of  the  hindered  rotor 

(7  9) 

functions  at  the  boundary  separating  reactants  and  products.  ’ This 

type  of  decoupling  in  NCC  has  been  shown  to  be  remarkably  effective  for 

(4  7) 

H + H,,  where  comparison  with  accurate  close-coupling  is  possible.  ’ 

This  success  is  partly  due  to  the  use  of  a switching  angle  to  reorient  the 
Z-axis  from  reactants  to  products  during  the  reaction. 


B.  The  Representation  of  Potential  Energy 

Surfaces  in  Natural  Collision  Coordinates 


The  convenience  of  NCC  in  studying  reactive  systems  arises  from 
the  use  of  a single  variable,  s,  to  follow  the  system  from  reactants  to 
products,  and  internal  (vibration-rotation)  coordinates  p,  > that  vary 
smoothly  from  reactants  to  products,  s is  taken  to  be  the  position  on  a 
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reference  curve  defined  in  the  eollinear  plane  corresponding  to  an  instan- 
taneous position  of  the  system  point  representing  the  reacting  molecules. 

It  is  through  the  potential  energy  decomposition 

V(s,p,y)  - V^s)  + V.,(p;s)  + V j (y  ; s)  + V^(p,y;s)  (15) 

that  the  concept  ot  a reaction  path  is  introduced.  is  defined  by  specifying 

P “ r -r  and  v in  a constant  s-plane.  By  connecting  the  values  of  V,  ob- 
o c o o 1 

tained  in  this  manner  and  plotting  them  as  a function  of  s we  obtain  a "reaction 
path  potential"  connecting  reactants  and  products,  i.e.,  a potential  energy 
correlation  diagram  for  the  adiabatic  electronic  state  of  interest.  Since  (p,y) 
are  orthogonal  coordinates,  once  is  specified,  V.,  and  then  represent 
stretching  and  bending  potentials  for  displacements  away  from  the  NCC  equi- 
librium geometry  (s«  p , v ) . V,  is  required  to  make  the  decomposition 
exact,  and  represents  vibration-rotation  interaction  in  the  potential  energy, 
which  for  some  systems  is  small  in  NCC.  This  representation  of  the  potential 
surface  is  very  convenient  because  the  translational,  vibrational  and  rotational 
potentials  can  often  be  represented  by  simple  analytic  functions.  This  simpli- 
fies studies  ot  the  effects  of  potential  surface  variation  on  the  dynamics  of 
a chemical  reaction.  The  separate  analytic  functions  defining  vibration  and 
rotation  represent  collective  motions  of  the  three-atom  system,  similar  to 
normal  coordinates.  For  H + H.?  at  s = 0 the  motions  in  NCC  correspond  to 
normal  coordinates. 

It  is  convenient  to  choose  V°  and  V^1  to  be  potentials  such  that  the 
resulting  hamiltonians  for  internal  motion  are  easily  soluble.  This  is 
because  we  need  eigenfunctions  of  internal  motion  to  expand  the  scattering 
wavef unction.  By  choosing  these  potentials  to  combine  with  kinetic  energy 
terms  to  give  Morse  oscillator  and  hindered  asymmetric  top  hamiltonians,  we 
then  have  a decomposition 

V(s,p,y)  = V1(s)  + V°(p;s)  + V°(y;s)  + V°(p,YJs)  (16) 

where  is  now  given  by 

V4  = V4  + ^V2~V2)  + (V3'V3) 


(17) 


We  have  i omul  that  Morse  and  hindered  rotor  eigenfunctions  are  very  useful 
representations  to  use  in  calculations.  Other  choices  of  representation, 
such  as  harmonic  oscillators  and  free  rotors,  allow  rapid  calculation  of 
nutrix  elements,  but  result  in  a much  larger  interaction  matrix  in  the 
scattering  equations,  since  they  do  not  represent  the  vibrational  and 
rotational  motions  as  well.  They  have  advantages  in  some  applications,  how- 
ever, particularly  the  harmonic  oscillator  functions,  since  they  can  be  used 
to  aid  convergence  when  continuum  vibrational  states  play  a role  in  the 
scattering.  Since  harmonic  oscillator  functions  and  corresponding  matrix 
elements  are  much  less  time  consuming  to  compute  than  those  of  Morse 
oscillators,  convergence  tests  and  studies  of  step-size  variations  are  more 
efficiently  performed  witli  these  simpler  functions. 

The  quantities  necessary  to  define  the  NCC  representation  of  the 
potential  surface  are  the  path  potential  V^,  and  quantities  parametric  in  s 
used  to  define  the  model  potentials  for  internal  motion.  These  are  r and 
y ^ (obtained  when  is  determined),  , the  local  Morse  exponent,  and  a 
functional  representation  of  the  local  bending  potential.  In  Figures  3 to  5 

we  illustrate  the  variation  of  V. , r , y , and  a for  several  reactions. 

1 o 'o  m 

Cubic  splines  are  used  to  represent  the  curves.  There  are  two  particularly 
useful  forms  for  the  bending  potential.  The  simple  form 

Vj(y)  = jkl-cos  2> ) Qg) 

and  the  much  more  flexible  Legendre  expansion 


V (y)  = l C P (cos  y)  (19) 

3 n n 

n 

The  coefficients  ~ and  are  parametric  in  s,  and  vary  smoothly  from  reac- 
tants to  products.  We  have  employed  the  first  form  extensively  for  our 
F + calculations.  The  second  form  is  required  for  a realistic  represen- 
tation of  a system  with  a bent  intermediate.  In  Figure  6 we  show  a fit 
of  for  F + In  Figure  7 we  show  a Legendre  fit  to  the  local  bending 


I 
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potential  tor  H + taken  from  the  potential  surface  of  Gauss.  It  should 
be  noted  that  the  motion  of  the  three  atom  system  in  NCC  depends  on  the  choice 
of  floating  origin  and  switching  angle.  Thus  the  potential  fits  will  vary 
with  this  choice. 

C.  Computational  Aspects 

Apart  from  the  representation  of  the  potential  energy  surface, 
scattering  calculations  can  be  considered  in  three  separate  steps.  First, 
the  matrix  elements  appearing  in  the  coupled  equations  must  be  computed. 

Second,  the  equations  must  be  integrated  along  the  reaction  coordinate  to 
obtain  a global  solution.  Finally,  the  physical  solution  is  obtained  by  the 
application  of  appropriate  boundary  conditions . ^ 

For  the  calculations  presented  in  this  report,  we  have  neglected 
vibration-rotation  coupling  in  the  potential,  i.e.,  we  set  = 0.  This  is 
consistent  with  our  treatment  of  the  kinetic  energy,  and  allows  us  to  vary 
the  vibrational  and  rotational  parts  of  the  problem  independently.  It  should 
be  pointed  out  that  this  pertains  to  NCC  only;  the  internuclear  coordinates 
are  still  coupled  in  a cartesian  space. 

Since  we  use  a vibrational  basis  of  Morse  oscillators  all  vibrational 
matrix  elements  must  be  evaluated  numerically.  We  use  standard  Simpson  quad- 
rature formulas  for  this  purpose.  The  rotational  matrix  elements  are  deter- 
mined analytically  in  a basis  of  free  rotors,  and  the  hindered  rotor  func- 
tions are  then  obtained  by  diagonalization  using  at  present  a Jacobi  algorithm. 
The  time  required  to  obtain  both  vibrational  and  rotational  matrix  elements 
for  a system  like  H + 0^  is  about  1000  CP  seconds  on  a CDC  7600,  assuming  16 
vibrational  states  with  20  rotational  levels.  We  have  found  it  more  efficient 
to  discard  the  basis  functions  and  retain  only  the  matrix  elements  used  for 
scattering.  This  is  because  of  the  large  amount  of  space  required  to  store 
the  numerical  functions.  The  matrix  elements  are  computed  once,  and  used 
for  all  energies,  and  for  all  angular  momenta  if  centrifugal  decoupling  is 
employed. 

Two  versions  of  the  scattering  code,  REACTOR,  were  employed.  The 
first,  an  initial-value  method  using  basis  sets  that  do  not  vary  within  a 


small  region  (sector)  along  the  reaction  coordinate,  using  the  exponential 
sector  transformation  ot  Redmon  and  Wyatt.  This  method  had  been  previously 
applied  to  the  H + , H + D and  F + ll.;  reactions,  and  represents  the 

sector  overlap  matrix  required  for  the  transformation  as 

-i6P  /ft 

_ s -i6P  /h  (20) 

S = e = e s 

This  use  of  the  translation  operator  requires  that  we  compute  the  energy- 
independent  first  derivative  matrix  P ^ . It  should  be  noted  that  we  are 
making  an  approximation,  in  representing  a matrix  exponential  as  the  expo- 
nential of  a matrix,  that  becomes  exact  in  the  limit  of  a complete  basis. 
Since  P is  antisymmetric  and  traceless,  S is  orthogonal,  and  flux  is  con- 
served  across  sector  boundaries.  Furthermore,  using  the  Light-Walker  method 
of  wavefunction  scaling^^  with  the  sector  method  described  above  (i.e., 
without  explicitly  including  first  and  second  derivative  couplings)  results 
in  a symmetric  interaction  matrix.  One  would  expect  then  that  a unitary 
S-matrix  would  result,  and  this  is  what  we  observe.  In  most  instances  total 
probability  is  unity  to  seven  or  eight  digits. 

The  main  difficulty  with  this  procedure  is  the  assumption  of  com- 
pleteness used  in  the  exponential  matrix  evaluation.  Excellent  results  were 
obtained  for  H + H^,  but  F + H.(  presented  difficulties  in  trying  to  achieve 
convergence,  although  the  probabilities  themselves  appeared  reasonable.  We 

therefore  modified  REACTOR  by  using  the  Light-Walker  R-matrix  propagation 

(12) 

method.  All  matrix  element  evaluations  were  the  same,  except  that  in- 

stead of  using  the  exponential  transformation  method,  requiring  first- 
derivatives,  direct  numerical  overlaps  are  used.  This  generally  is  expected 
to  degrade  flux  conservation  when  an  incomplete  basis  is  used,  but  because 
of  the  way  in  which  the  transformation  is  constructed  in  the  R-matrix  method, 
i.e.,  by  replacing  matrix  inverses  by  transposes  (again,  this  is  exact  in  the 
limit  of  a complete  basis)  flux  is  satisfactorily  conserved.  In  our  F + 
calculations  using  the  R-matrix  method,  flux  was  conserved  to  a few  tenths 


n 


of  one  percent.  While  this  is  not  as  impressive  as  the  conservation  ob- 
tained with  the  exponential  sector  transformations,  it  is  satisfactory.  The 
R-matrix  method  appears  more  stable  with  respect  to  convergence  of  the  vibra- 
tional basis.  There  is  no  apparent  reason  why  the  exponential  method  and 
the  R-matrix  method  should  have  the  same  rate  of  convergence,  and  in  fact 
our  difficulty  with  the  exponential  transformation  may  stem  from  the  slow 
convergence  of  our  representation  of  the  matrix  of  the  translation  operator. 

It  is  apparently  much  slower  for  F + than  the  rate  of  convergence  of  the 
coupled-channel  expansion  of  the  wavef unction. 

The  initial-value  method  is  inherently  unstable  to  closed  channels. 
In  practice,  one  observes  that  total  probabilities  of  unity  can  be  obtained 
if  closed-channels  are  not  carried  to  asymptotic  limits  of  the  reaction  co- 
ordinate, but  are  picked  up  and  dropped  at  appropriate  places  during  the 
course  of  the  integration.  Experience  has  taught  us  how  to  do  this  effec- 
tively. The  R-matrix  method  does  not  have  this  difficulty,  and  closed 
channels  can  be  carried  to  asymptotic  regions.  We  feel  that  the  sensitivity 
of  the  initial  value  approach  to  closed-channels  can  be  viewed  in  a positive 
light,  since  when  conservation  is  achieved  one  can  perhaps  feel  confident 
about  the  stability  of  the  integration.  Because  closed  channels  are  not 
carried  to  asymptotic  limits,  our  current  implementation  of  the  initial 
value  method  is  somewhat  more  efficient  than  our  R-matrix  code.  However, 
there  is  no  reason  why  the  R-matrix  method  cannot  be  similarly  coded. 

Most  of  the  applications  presented  in  this  report  were  obtained 
using  the  R-matrix  version  of  REACTOR.  The  most  time  consuming  calculations 
involved  the  60  channel  runs  for  F + , carried  out  on  the  CDC  6600  at  the 

University  of  Texas  at  Austin  in  collaboration  with  Professor  Wyatt.  We 
used  280  steps  over  the  reaction  coordinate,  and  typical  execution  times 
were  1500  CP  seconds  per  E and  J.  We  could  have  saved  time  by  using  200 
steps,  with  only  slight  loss  in  accuracy.  Most  of  the  time  is  spent  in 

matrix  multiplication  and  diagonalization  routines.  A Jacobi  method  was 
used  for  diagonalization,  and  is  known  to  be  inefficient  for  large  matrices. 
Code  development  currently  in  progress  (including  optimum  step-size  selection) 
should  enable  us  to  obtain  similar  results  in  about  600-800  CP  seconds.  This 


corresponds  to  about  100  CP  seconds  per  E,  .1  on  a CDC  7600.  This  means  that 
the  information  required  to  obtain  reliable  state-to-state  rate  constants 
tor  reactions  like  F + II.,,  0 + H,  and  H + 0 can  be  obtained  for  a few  hours 
of  computation  on  a CDC  7600. 


D.  Recent  Results  tor  the  F + U2  Reaction 

The  R-matrix  version  of  REACTOR  has  been  applied  to  the  H + H,,  re- 
actions, and  reaction  probabilities  are  in  good  agreement  witli  results  ob- 
(4  5) 

talned  previously.  * An  extensive  set  of  calculations  wore  undertaken  to 
study  the  effects  ot  varying  the  bending  potential,  V^,  on  computed  reaction 
probabilities.  This  type  of  Investigation  is  expected  to  yield  information 
concerning  the  importance  of  nonlinear  configurations  in  determining  the 
product  energy  disposal  from  reactive  encounters.  Similar  collinear  quantum 
calculations  investigated  the  effects  of  variation  of  the  F + H~  vibrational 
potential  on  reaction  probabilities.^^  It  was  found  that  the  presence  of 
wells  in  the  adiabatic  vibrational  correlation  diagram  were  responsible  for 
the  low-energy  population  inversion  observed  in  the  reactive  probability 
curves.  These  wells  are  associated  with  the  local  expansion  of  the  collinear 
potential  surface  in  the  downhill  part  of  the  reaction  path.  This  is  re- 
flected in  a drop  in  the  value  of  the  local  Morse  exponent  used  to  define  the 
vibrational  potential,  and  Indicates  that  an  accurate  knowledge  of  the  features 
of  the  potential  surface  in  this  region  is  crucial  for  obtaining  state-to-state 
reaction  probabilities. 

A model  potential  of  the  form  given  by  Equation  (18)  was  used  for 

Vo 

the  bending  potential  in  these  investigations.  We  tirst  represented  y by 
the  simple  parametric  form 


Vo  _ A 

2 = -BX.  BX  (21) 

e +e 


where  X - |s-s  I.  By  varying  the  parameters  in  the  expression,  the 
1 max  1 

hindered  rotor  eigenvalues  obtained  with  the  potential  of  Equation  (18)  can 

be  made  to  have  a maximum  value  at  s = s . A is  chosen  so  that  the  lowest 

max 


(j=0)  eigenvalue  has  a preselected  maximum.  We  initially  chose  this  value 

to  correspond  to  the  zero-point  bending  energy  of  the  Muckerman  II  F + H.  sur- 

(14)  1 

face.  Two  values  were  used  for  s , s = 0.49  corresponding  to  the  F + H„ 

max  2 

saddle  point  on  the  Muckerman  V surface,  and  s = 0 corresponding  to  a 

point  on  the  downhill  part  of  the  reaction  path.  The  latter  value  was  chosen 
because  by  placing  the  maximum  in  the  bending  zero-point  energy  at  a value 
of  s where  both  the  potential  and  the  vibrational  energy  are  dropping 
rapidly  (i.e.,  we  are  going  into  a vibrational  well)  we  expect  to  minimize 
the  effect  of  the  bending  energy  on  the  reaction  dynamics.  A final  modifi- 
cation of  the  bending  potential,  now  being  used  in  extensive  calculations, 

Vo 

is  to  replace  equation  21  with  a locally  determined  fit  to  — which  repro- 
duces the  j = 0 eigenvalue  for  the  Muckerman  V surface  at  many  values  of  s. 
This  should  allow  us  to  model  the  dynamics  on  this  surface  rather  well.  It 
should  be  mentioned  that  the  reason  for  not  using  an  accurate  Legendre  ex- 
pansion such  as  Equation  (19)  is  our  desire  to  eliminate  the  complication 
of  the  term  in  the  potential,  and  a of  the  form  of  Equation  (19)  results 
in  a bending  potential  that  is  much  too  hard,  requiring  that  be  included 
in  order  that  the  rotational  eigenvalues  should  not  be  too  large. 

The  extensive  results  we  are  getting  for  this  system  are  to  be  pub- 
lished in  the  scientific  literature , ^^)  aruj  therefore  will  only  be  briefly 
discussed  here.  Our  calculations  for  this  particular  study  involve  10 
vibrational  states,  and  up  to  12  rotational  states  in  the  4 open  product 
vibrational  levels,  for  a total  of  60  coupled  channels.  Previous  work  has 
shown  that  this  basis  is  sufficiently  converged  for  our  purpose.  In  Figure  8 
we  attempt  to  represent  in  perspective  the  information  that  is  directly  ob- 
tained from  calculations  of  this  type.  Scattering  probabilities  are  obtained 
as  a function  of  total  angular  momentum  J and  total  energy  E,  from  which  it 
is  possible  to  obtain  cross-sections  and  finally  rate  constants.  We  show 
in  this  figure  the  probability  for  ground  state  reacting  with  F atoms  to 
produce  vibrationally  excited  HF  in  V'  = 2.  This  figure  corresponds  to 

s = 0 in  Equation  (21).  The  low  energy  probability  curve  for  J * 0 (vs.  E) 
max 

shows  a pronounced  maximum  around  .34  eV.  Going  out  in  total  J (at  this 
energy)  the  probability  monotonically  decreases.  If  we  choose  a higher  energy 
and  scan  out  in  total  J,  we  notice  instead  that  the  probability  rises  to  a 


Ji 


maximum  at  a nonzero  value  of  J before  exhibiting  monotonical ly  decreasing 
behavior.  The  explanation  for  this  lies  in  the  shape  of  the  J = 0 probability 
vs  energy.  The  sharp  maximum  at  .34  eV  is  the  result  of  a Feshbach  resonance 
process^*^  involving  quasibound  states  associated  with  wells  in  the  dis- 
tortion potentials  mentioned  above.  For  -I%0,  the  probability  curves  (vs  E) 
are  similar,  but  the  onset  ot  resonance  behavior  occurs  at  higher  total  energy 
than  at  J = 0 due  to  the  shift  of  reaction  threshold.  These  are  quantum 
effects,  and  total  cross-sections  from  our  previous  calculations  indicate 
that  they  are  not  completely  removed  by  averaging. ^ ^ 

In  Figure  9 the  probabilities  are  shown  in  more  detail  for  two 
energies.  In  Figure  10  the  probabilities  are  plotted  for  .1=0  and  J = 10 
as  a function  of  energy.  Also  shown  in  these  figures  are  probabilities  for 
reaction  into  the  V'  = 3 manifold  of  HF. 

Calculations  in  which  s = 0.49  result  in  very  small  reaction 

probabilities  for  J = 0.  One  maximum  of  .07  occurs  at  .35  eV,  with  another 

of  .08  at  .43  eV.  The  reaction  probabilities  are  thus  drastically  altered 

Vo 

by  moving  the  maximum  in  — along  the  reaction  path.  With  this  potential, 
we  find  a large  reaction  probability  only  for  J>0,  where  we  again  observe  the 
resonance  behavior  in  the  probabilities. 

(4) 

The  first  3-1)  quantum  calculations  that  we  reported  for  F + H., 
were  made  with  a very  limited  rotor  basis  (2  rotors  per  vibrational  level) 
and  an  integration  method  that  explicitly  included  first  and  second  deriva- 
tives of  the  basis  functions.  These  results  showed  the  population  inversion 
into  V'  = 2 of  HF,  and  subsequent  work  lias  verified  the  essential  correctness 
of  the  probability  curves  reported  there.  We  have  performed  similar  calcu- 
lations using  the  present  model  bending  potential,  with  s = 0.  The  results 
are  shown  in  Figure  11,  and  are  not  significantly  altered  from  the  results 
of  Ref.  4,  even  though  the  bending  potential  in  that  case  was  an  expansion 
of  the  form 


V3(y)  - ~Vo(s)(l-cos2Y)  + A(s)£  C„Pn (cosy)+B(s)  (2i 

n=0 

in  which  the  I.egendre  coefficients  were  fixed  at  one  value  of  s (to  set  the 
shape  of  the  potential)  with  Vo,  A and  1)  determined  locally.  These  studies 


indicate  ttiat  qualitatively  correct  reaction  probabilities  for  vibrational 

processes  may  be  obtained  for  some  systems  with  severely  truncated  rotor  bases. 

Vo 

Currently  we  are  using  a V ^ given  by  Equation  (18)  with  — as  shown 
in  Figure  b in  an  extensive  study  of  the  F + reaction,  in  which  state-to- 
state  reaction  rates  will  be  the  end  result. 


E.  Reaction  Probabilities  tor  the  H + 02  Reaction 


We  have  applied  the  code  REACTOR  to  the  H + 0.,  reaction  using  a 


fit  to  the  modified  LEPS  surface  ot  Gauss. 


A vibrational  basis  of  up 


to  16  Morse  oscillators  was  used  with  a limited  rotational  basis  of  only  2 

(4) 

hindered  rotors,  in  the  spirit  of  our  early  F + H.,  calculations.  Pre- 

liminary results  for  vibrationally  excited  0_  reacting  to  produce  ground 
state  OH  indicate  that  vibrational  energy  is  very  efficient  in  promoting  the 
reaction.  This  is  in  agreement  with  the  classical  trajectory  results  of 
Gauss.  We  have  attempted  to  determine  the  convergence  of  the  coupled 

channel  expansion,  but  we  are  presently  limited  by  code  dimensions  to  about 
70  channels,  and  the  behavior  observed  so  far  indicates  we  will  need  over 
100.  Flux  conservation  is  usually  better  than  l*  for  the  calculations  re- 
ported here.  When  we  tried  to  increase  the  rotational  basis,  flux  conserva- 
tion deteriorated  rapidly,  due  to  mixing,  and  probably  will  not  improve  until 
many  more  states  are  used.  A modified  version  of  the  program  employing  ex- 
tended core  storage  will  be  used  for  calculations  with  a larger  basis.  This 
will  allow  us  to  Investigate  the  role  of  reagent  rotational  energy  on  the 
scattering. 

In  Figure  12  we  present  reaction  probabilities  for  H + 0o(v)  at 
1.0  eV  total  energy.  It  is  evident  that  reagent  vibrational  energy  is  very 
effective  in  promoting  reaction.  Comparisons  of  these  figures  with  similar 
results  for  F + H9  in  Figure  9 show  significant  differences  in  the  dynamics 
between  the  two  reactions.  The  H + 0.,  probability  curves  have  more  struc- 
ture, probably  indicative  ot  complexing  associated  with  the  deep  well  in  the 
surface.  In  contrast,  the  F + H2  curves  are  relatively  simple. 

In  Table  II  we  present  some  probabilities  for  H + 02(v=4,  j=0)  at 
several  energies  and  total  angular  momenta.  The  structure  indicated  by  the 
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diversity  in  magnitude  of  these  results  illustrates  the  complexity  we  gene- 
rally find  in  the  energy  variation  of  the  scattering  probabilities.  By  con- 
trast, the  F + probabilities  are  quite  smooth.  Because  of  this  structure, 
it  will  be  necessary  to  solve  the  coupled  equations  over  a rather  fine 
energy  grid. 

We  find  that  reaction  from  v=0  in  0,,  is  small,  but  not  negligible. 
Gauss  found  essentially  no  classical  trajectories  that  were  reactive  when 
starting  out  in  v=0.  It  is  likely  that  some  of  ills  trajectories  that  were 
labeled  as  complex  would  actually  have  led  to  reaction  if  tie  could  have 
followed  them  to  completion,  lie  was  not  able  to  do  this  because  of  numerical 
difficulties  encountered  in  trying  to  integrate  such  trajectories. 

Another  interesting  feature  of  the  reaction  is  the  difference  between 
the  v=4  curve  and  the  curves  for  v=2  and  v=3  in  Figure  12.  The  v=4  curve 
grows  gradually  and  has  a very  broad  maximum,  while  the  others  drop  to  a low 
value  at  relatively  small  J.  This  indicates  that  the  large  cross-sections 
found  for  v>4  may  be  associated  witli  dynamical  effects  that  are  quite  different 
than  occur  when  less  vibrational  energy  is  available.  Gauss  ran  only  a few 
trajectories  tor  v<4  and  these  gave  no  reaction.  It  would  be  useful  to  see 
if  trajectories  for  v<4  might  indicate  a more  direct  mechanism,  with  less  com- 
plexing. 


F.  Gross-Sections  for  the  H + O2  Reaction 


The  probabilities  discussed  in  the  preceding  section  have  been  used 

to  estimate  total  cross-sections  for  h + 0.  at  1.0  eV  total  energy.  These  are 

^ o 2 

presented  in  Table  III.  Gauss  reported  a value  of  about  .35  A for  the  v=4 
vibrational  state  of  0^.  He  estimates  that  his  value  could  be  as  much  as  a 
factor  of  10  too  small,  a result  of  his  inability  to  follow  complex  trajec- 
tories. Our  quantum  cross-section  is  about  10  times  larger  than  his,  which 
is  consistent  with  his  estimate  of  the  number  of  complex  trajectories  that 
might  lead  to  reaction.  It  is  necessary  for  us  to  obtain  cross-sections  at 
other  energies  before  a rate-constant  can  be  determined  for  this  reaction. 

The  rate-constants  obtained  by  Gauss  are  lower  than  experiment,  and  the  larger 
cross-sections  that  seem  to  be  resulting  from  the  quantum  scattering  are 
evidently  going  to  Improve  this,  although  the  threshold  behavior  will  have  to 
be  carefully  investigated. 


1 — 
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In  Table  IV  we  present  some  inelastic  (nonreactive)  partial  cross- 
sections  QJ  which  indicate  that  inelastic  H + 0^  encounters  tend  to  produce 
vibrationally  excited  which  in  turn  can  react  rapidly  with  hydrogen. 

This  could  complicate  the  analysis  of  experiments  in  which  state-selection 
is  not  employed. 

The  LEPS  potential  surface  used  for  these  calculations  does  not 

produce  the  known  force  constants  for  HC^.  and  since  these  calculations  were 

begun.  Gauss  has  modified  the  surface  to  reproduce  these  force  constants  to 
(18) 

about  20%.  ' In  addition,  an  extensive  ab-initio  surface  obtained  with 

(19) 

large-scale  Cl  calculations  has  become  available.  This  latter  surface 

should  be  employed  in  any  extensive  studies  of  the  H + O2  reaction. 
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Table  1 - Combustion  Reactions  of  Interest 


H + 0 -►  OH  + 0 

0 + H2  -*■  OH  + H 
0 + N2  -*•  ON  + N 
N + 02  -*■  NO  + 0 
B + OH  •+  BO  + H 
B + 02  -*•  BO  + 0 

c + o9-*-co  + o 

C + OH  ->  CO  + H 
CN  + 0 ->  CO  + N 
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Table  II  - Reaction  Probabilities  for 

H + 02(v=4,  j =0)  -+  0 + OH(v=0,  j=all) 


etotal  (eV) 

J=0 

5 

10 

15 

20 

.85 

.161 

.143 

.080 

0. 

0. 

.90 

.002 

.236 

.023 

.093 

.020 

.95 

.097 

.013 

.491 

.040 

.007 

1.00 

.113 

.249 

.406 

.083 

.004 

1.05 

.090 

1.10 

.146 

Table  IV  - Nonreactive  partial  cross-sections* 


for  H + 

02(v=0;  j=0)  - 

H + °2(V»J 

V 

. 85eV , J=5 

. 9eV , J=10 

l.OeV,  : 

0 

.01 

.11 

.09 

1 

.15 

.15 

.17 

2 

.08 

.06 

.12 

3 

.05 

.24 

.38 

4 

.04 

.03 

.09 

* Values  are  in  A . 
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CAPTIONS 

1 The  natural  coordinate  representation  of  the  collinear  plane. 
The  reaction  coordinate  s is  measured  along  the  reference 
curve  RC,  with  s=0  corresponding  to  its  intersection  with  SS. 
Also  shown  is  a planar  reaction  path. 

2 Deviation  from  linearity  in  a plane  of  constant  s.  The  region 

between  v and  the  m-axis  must  be  avoided  to  avoid  double- 
m 

counting  of  spatial  configurations  of  the  3-atom  system. 

3 Reaction  Path  potentials  for  H + H2>  F + H2>  and  H + 02 
obtained  from  empirical  LEPS  surfaces. 

4 Exponential  parameter  u used  to  define  local  Morse  potentials 
for  F + and  h + 09. 

5 The  bending  angle  y o and  scaled  NCC  equilibrium  distance  rQ 
which  define  the  reaction  path  for  the  H + C>2  reaction. 

Vo 

6 A local  fit  to  — used  to  parametrize  the  bending  potential 
for  F + H2- 

7 Ten-term  Legendre  fits  to  the  local  bending  potential  for 
H + C>2 » shown  for  three  values  of  s. 

8 A perspective  representation  of  the  P reaction  probability 
for  F + H2  as  a function  of  both  E and  J. 

F + H2  reaction  probabilities  as  a function  of  J at  two 
energies. 


Figure  9 
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FIGURE  CAPTIONS 


Figure  1 


Figure  2 


Figure  3 


Figure  4 


Figure  5 


Figure  6 


Figure  7 


The  natural  coordinate  representation  of  the  collinear  plane. 
The  reaction  coordinate  s is  measured  along  the  reference 
curve  RC,  with  s=0  corresponding  to  its  intersection  with  SS. 
Also  shown  is  a planar  reaction  path. 


Deviation  from  linearity  in  a plane  of  constant  s.  The  legion 

between  v and  the  m-axis  must  be  avoided  to  avoid  double- 
m 

counting  of  spatial  configurations  of  the  3-atom  system. 


Reaction  Path  potentials  for  H + F + I^,  and  H + 

obtained  from  empirical  LEPS  surfaces. 

Exponential  parameter  a used  to  define  local  Morse  potentials 
for  F + and  H + C^. 

The  bending  angle  y^  and  scaled  NCC  equilibrium  distance  r^ 
which  define  the  reaction  path  for  the  h + reaction. 


A local  fit  to  — used  to  parametrize  the  bending  potential 
for  F + 


Ten-term  Legendre  fits  to  the  local  bending  potential  for 
H + 0^,  shown  for  three  values  of  s. 


Figure  8 A perspective  representation  of  the  P ^ reaction  probability 

for  F + as  a function  of  both  E and  J. 

Figure  9 F + H„  reaction  probabilities  as  a function  of  J at  two 

energies . 


Figure  10  F + reaction  probabilities  as  a function  of  E for 

.1=0  and  .1  = 10. 

Figure  11  F + reaction  probabilities  as  a function  of  E for  J-0, 

obtained  with  a truncated  rotational  basis. 

Figure  12  H + 0^  reaction  probabilities  as  a function  of  J for  E-1.0  eV, 

obtained  with  a truncated  rotational  basis. 
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